Catalytic Metal Ion-Substrate Coordination during Nonenzymatic RNA Primer Extension.

Nonenzymatic template-directed RNA copying requires catalysis by divalent metal ions. The primer extension reaction involves the attack of the primer 3'-hydroxyl on the adjacent phosphate of a 5'-5'-imidazolium-bridged dinucleotide substrate. However, the nature of the interaction of the catalytic metal ion with the reaction center remains unclear. To explore the coordination of the catalytic metal ion with the imidazolium-bridged dinucleotide substrate, we examined catalysis by oxophilic and thiophilic metal ions with both diastereomers of phosphorothioate-modified substrates. We show that Mg2+ and Cd2+ exhibit opposite preferences for the two phosphorothioate substrate diastereomers, indicating a stereospecific interaction of the divalent cation with one of the nonbridging phosphorus substituents. High-resolution X-ray crystal structures of the products of primer extension with phosphorothioate substrates reveal the absolute stereochemistry of this interaction and indicate that catalysis by Mg2+ involves inner-sphere coordination with the nonbridging phosphate oxygen in the pro-SP position, while thiophilic cadmium ions interact with sulfur in the same position, as in one of the two phosphorothioate substrates. These results collectively suggest that during nonenzymatic RNA primer extension with a 5'-5'-imidazolium-bridged dinucleotide substrate the interaction of the catalytic Mg2+ ion with the pro-SP oxygen of the reactive phosphate plays a crucial role in the metal-catalyzed SN2(P) reaction.

From Enigma to Revelation: Unravelling Biological Functions of Ubiquitous Small Ribozymes.

RNA, widely recognized as an information-carrier molecule, is capable of catalyzing essential biological processes through ribozymes. Despite their ubiquity, specific functions in a biological context and phenotypes based on the ribozymes' activity are often unknown. Here, we present the discovery of a subgroup of minimal HDV-like ribozymes, which reside 3' to viral tRNAs and appear to cleave the 3'-trailers of viral premature tRNA transcripts. This proposed tRNA-processing function is unprecedented for any ribozymes, thus, we designate this subgroup as theta ribozymes. Most theta ribozymes were identified in Caudoviricetes bacteriophages, the main constituent (>90%) of the mammalian gut virome. Intriguingly, our findings further suggest the involvement of theta ribozymes in the transition of certain bacteriophages between distinct genetic codes, thus possibly contributing to the phage lysis trigger. Our discovery expands the limited repertoire of biological functions attributed to HDV-like ribozymes and provides insights into the fascinating world of RNA catalysis.

High-Fidelity RNA Copying via 2',3'-Cyclic Phosphate Ligation.

Templated ligation offers an efficient approach to replicate long strands in an RNA world. The 2',3'-cyclic phosphate (>P) is a prebiotically available activation that also forms during RNA hydrolysis. Using gel electrophoresis and high-performance liquid chromatography, we found that the templated ligation of RNA with >P proceeds in simple low-salt aqueous solutions with 1 mM MgCl2 under alkaline pH ranging from 9 to 11 and temperatures from -20 to 25 °C. No additional catalysts were required. In contrast to previous reports, we found an increase in the number of canonical linkages to 50%. The reaction proceeds in a sequence-specific manner, with an experimentally determined ligation fidelity of 82% at the 3' end and 91% at the 5' end of the ligation site. With splinted oligomers, five ligations created a 96-mer strand, demonstrating a pathway for the ribozyme assembly. Due to the low salt requirements, the ligation conditions will be compatible with strand separation. Templated ligation mediated by 2',3'-cyclic phosphate in alkaline conditions therefore offers a performant replication and elongation reaction for RNA on early Earth.

Multiple structural flavors of RNase P in precursor tRNA processing.

The precursor transfer RNAs (pre-tRNAs) require extensive processing to generate mature tRNAs possessing proper fold, structural stability, and functionality required to sustain cellular viability. The road to tRNA maturation follows an ordered process: 5'-processing, 3'-processing, modifications at specific sites, if any, and 3'-CCA addition before aminoacylation and recruitment to the cellular protein synthesis machinery. Ribonuclease P (RNase P) is a universally conserved endonuclease in all domains of life, performing the hydrolysis of pre-tRNA sequences at the 5' end by the removal of phosphodiester linkages between nucleotides at position -1 and +1. Except for an archaeal species: Nanoarchaeum equitans where tRNAs are transcribed from leaderless-position +1, RNase P is indispensable for life and displays fundamental variations in terms of enzyme subunit composition, mechanism of substrate recognition and active site architecture, utilizing in all cases a two metal ion-mediated conserved catalytic reaction. While the canonical RNA-based ribonucleoprotein RNase P has been well-known to occur in bacteria, archaea, and eukaryotes, the occurrence of RNA-free protein-only RNase P in eukaryotes and RNA-free homologs of Aquifex RNase P in prokaryotes has been discovered more recently. This review aims to provide a comprehensive overview of structural diversity displayed by various RNA-based and RNA-free RNase P holoenzymes towards harnessing critical RNA-protein and protein-protein interactions in achieving conserved pre-tRNA processing functionality. Furthermore, alternate roles and functional interchangeability of RNase P are discussed in the context of its employability in several clinical and biotechnological applications. This article is categorized under: RNA Processing > tRNA Processing RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

New Insights into the Dependence of CPEB3 Ribozyme Cleavage on Mn2+ and Mg2.

CPEB3 ribozyme is a self-cleaving RNA that occurs naturally in mammals and requires divalent metal ions for efficient activity. Ribozymes exhibit preferences for specific metal ions, but the exact differences in the catalytic mechanisms of various metal ions on the CPEB3 ribozyme remain unclear. Our findings reveal that Mn2+ functions as a more effective cofactor for CPEB3 ribozyme catalysis compared to Mg2+, as confirmed by its stronger binding affinity to CPEB3 by EPR. Cleavage assays of CPEB3 mutants and molecular docking analyses further showed that excessive Mn2+ ions can bind to a second binding site near the catalytic site, hindering CPEB3 catalytic efficiency and contributing to the Mn2+ bell-shaped curve. These results implicate a pivotal role for the local nucleobase-Mn2+ interactions in facilitating RNA folding and modulating the directed attack of nucleophilic reagents. Our study provides new insights and experimental evidence for exploring the divalent cation dependent cleavage mechanism of the CPEB3 ribozyme.

Hydrolytic endonucleolytic ribozyme (HYER) is programmable for sequence-specific DNA cleavage.

Ribozymes are catalytic RNAs with diverse functions including self-splicing and polymerization. This work aims to discover natural ribozymes that behave as hydrolytic and sequence-specific DNA endonucleases, which could be repurposed as DNA manipulation tools. Focused on bacterial group II-C introns, we found that many systems without intron-encoded protein propagate multiple copies in their resident genomes. These introns, named HYdrolytic Endonucleolytic Ribozymes (HYERs), cleaved RNA, single-stranded DNA, bubbled double-stranded DNA (dsDNA), and plasmids in vitro. HYER1 generated dsDNA breaks in the mammalian genome. Cryo-electron microscopy analysis revealed a homodimer structure for HYER1, where each monomer contains a Mg2+-dependent hydrolysis pocket and captures DNA complementary to the target recognition site (TRS). Rational designs including TRS extension, recruiting sequence insertion, and heterodimerization yielded engineered HYERs showing improved specificity and flexibility for DNA manipulation.

Chemically Cross-Linked Hammerhead Ribozyme as an Efficient RNA Interference Tool.

RNA-cleaving ribozymes are promising candidates as general tools of RNA interference (RNAi) in gene manipulation. However, compared with other RNA systems, such as siRNA and CRISPR technologies, the ribozyme tools are still far from broad applications on RNAi due to their poor performance in the cellular context. In this work, we report an efficient RNAi tool based on chemically modified hammerhead ribozyme (HHR). By the introduction of an intramolecular linkage into the minimal HHR to reconstruct the distal interaction within the tertiary ribozyme structure, this cross-linked HHR exhibits efficient RNA substrate cleavage activities with almost no sequence constraint. Cellular experiments suggest that both exogenous and endogenous RNA expression can be dramatically knocked down by this HHR tool with levels comparable to those of siRNA. Unlike the widely applied protein-recruiting RNA systems (siRNA and CRISPR), this ribozyme tool functions solely on RNA itself with great simplicity, which may provide a new approach for gene manipulation in both fundamental and translational studies.

Deep generative design of RNA family sequences.

RNA engineering has immense potential to drive innovation in biotechnology and medicine. Despite its importance, a versatile platform for the automated design of functional RNA is still lacking. Here, we propose RNA family sequence generator (RfamGen), a deep generative model that designs RNA family sequences in a data-efficient manner by explicitly incorporating alignment and consensus secondary structure information. RfamGen can generate novel and functional RNA family sequences by sampling points from a semantically rich and continuous representation. We have experimentally demonstrated the versatility of RfamGen using diverse RNA families. Furthermore, we confirmed the high success rate of RfamGen in designing functional ribozymes through a quantitative massively parallel assay. Notably, RfamGen successfully generates artificial sequences with higher activity than natural sequences. Overall, RfamGen significantly improves our ability to design functional RNA and opens up new potential for generative RNA engineering in synthetic biology.

An RNA Motif That Enables Optozyme Control and Light-Dependent Gene Expression in Bacteria and Mammalian Cells.

The regulation of gene expression by light enables the versatile, spatiotemporal manipulation of biological function in bacterial and mammalian cells. Optoribogenetics extends this principle by molecular RNA devices acting on the RNA level whose functions are controlled by the photoinduced interaction of a light-oxygen-voltage photoreceptor with cognate RNA aptamers. Here light-responsive ribozymes, denoted optozymes, which undergo light-dependent self-cleavage and thereby control gene expression are described. This approach transcends existing aptamer-ribozyme chimera strategies that predominantly rely on aptamers binding to small molecules. The optozyme method thus stands to enable the graded, non-invasive, and spatiotemporally resolved control of gene expression. Optozymes are found efficient in bacteria and mammalian cells and usher in hitherto inaccessible optoribogenetic modalities with broad applicability in synthetic and systems biology.

Unusual Base Pair between Two 2-Thiouridines and Its Implication for Nonenzymatic RNA Copying.

2-Thiouridine (s2U) is a nucleobase modification that confers enhanced efficiency and fidelity both on modern tRNA codon translation and on nonenzymatic and ribozyme-catalyzed RNA copying. We have discovered an unusual base pair between two 2-thiouridines that stabilizes an RNA duplex to a degree that is comparable to that of a native A:U base pair. High-resolution crystal structures indicate similar base-pairing geometry and stacking interactions in duplexes containing s2U:s2U compared to those with U:U pairs. Notably, the C═O···H-N hydrogen bond in the U:U pair is replaced with a C═S···H-N hydrogen bond in the s2U:s2U base pair. The thermodynamic stability of the s2U:s2U base pair suggested that this self-pairing might lead to an increased error frequency during nonenzymatic RNA copying. However, competition experiments show that s2U:s2U base-pairing induces only a low level of misincorporation during nonenzymatic RNA template copying because the correct A:s2U base pair outcompetes the slightly weaker s2U:s2U base pair. In addition, even if an s2U is incorrectly incorporated, the addition of the next base is greatly hindered. This strong stalling effect would further increase the effective fidelity of nonenzymatic RNA copying with s2U. Our findings suggest that s2U may enhance the rate and extent of nonenzymatic copying with only a minimal cost in fidelity.

Co-transcriptional folding of the glmS ribozyme enables a rapid response to metabolite.

The glmS ribozyme riboswitch, located in the 5' untranslated region of the Bacillus subtilis glmS messenger RNA (mRNA), regulates cell wall biosynthesis through ligand-induced self-cleavage and decay of the glmS mRNA. Although self-cleavage of the refolded glmS ribozyme has been studied extensively, it is not known how early the ribozyme folds and self-cleaves during transcription. Here, we combine single-molecule fluorescence with kinetic modeling to show that self-cleavage can occur during transcription before the ribozyme is fully synthesized. Moreover, co-transcriptional folding of the RNA at a physiological elongation rate allows the ribozyme catalytic core to react without the downstream peripheral stability domain. Dimethyl sulfate footprinting further revealed how slow sequential folding favors formation of the native core structure through fraying of misfolded helices and nucleation of a native pseudoknot. Ribozyme self-cleavage at an early stage of transcription may benefit glmS regulation in B. subtilis, as it exposes the mRNA to exoribonuclease before translation of the open reading frame can begin. Our results emphasize the importance of co-transcriptional folding of RNA tertiary structure for cis-regulation of mRNA stability.

A comparative survey of the influence of small self-cleaving ribozymes on gene expression in human cell culture.

Self-cleaving ribozymes are versatile tools for synthetic biologists when it comes to controlling gene expression. Up to date, 12 different classes are known, and over the past decades more and more details about their structure, cleavage mechanisms and natural environments have been uncovered. However, when these motifs are applied to mammalian gene expression constructs, the outcome can often be unexpected. A variety of factors, such as surrounding sequences and positioning of the ribozyme influences the activity and hence performance of catalytic RNAs. While some information about the efficiency of individual ribozymes (each tested in specific contexts) is known, general trends obtained from standardized, comparable experiments are lacking, complicating decisions such as which ribozyme to choose and where to insert it into the target mRNA. In many cases, application-specific optimization is required, which can be very laborious. Here, we systematically compared different classes of ribozymes within the 3'-UTR of a given reporter gene. We then examined position-dependent effects of the best-performing ribozymes. Moreover, we tested additional variants of already widely used hammerhead ribozymes originating from various organisms. We were able to identify functional structures suited for aptazyme design and generated highly efficient hammerhead ribozyme variants originating from the human genome. The present dataset will aide decisions about how to apply ribozymes for affecting gene expression as well as for developing ribozyme-based switches for controlling gene expression in human cells.

Phosphonate and Thiasugar Analogues of Glucosamine-6-phosphate: Activation of the glmS Riboswitch and Antibiotic Activity.

The glmS riboswitch is a motif found in 5'-untranslated regions of bacterial mRNA that controls the synthesis of glucosamine-6-phosphate (GlcN6P), an essential building block for the bacterial cell wall, by a feedback mechanism. Activation of the glmS riboswitch by GlcN6P mimics interferes with the ability of bacteria to synthesize its cell wall. Accordingly, GlcN6P mimics acting as glmS activators are promising candidates for future antibiotic drugs that may overcome emerging bacterial resistance against established antibiotics. We describe the synthesis of a series of phosphonate mimics of GlcN6P as well as the thiasugar analogue of GlcN6P. The phosphonate mimics differ in their pKa value to answer the question of whether derivatives with a pKa matching that of GlcN6P would be efficient glmS activators. We found that all derivatives activate the riboswitch, however, less efficiently than GlcN6P. This observation can be explained by the missing hydrogen bonds in the case of phosphonates and is valuable information for the design of future GlcN6P mimics. The thiasugar analogue of GlcN6P on the other hand turned out to be a glmS riboswitch activator with the same activity as the natural metabolite GlcN6P. The nonphosphorylated thiasugar displayed antimicrobial activity against certain bacilli. Therefore, the compound is a promising lead structure for the development of future antibiotics with a potentially novel mode of action.

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time.

Folding of ribozymes into well-defined tertiary structures usually requires divalent cations. How Mg2+ ions direct the folding kinetics has been a long-standing unsolved problem because experiments cannot detect the positions and dynamics of ions. To address this problem, we used molecular simulations to dissect the folding kinetics of the Azoarcus ribozyme by monitoring the path each molecule takes to reach the folded state. We quantitatively establish that Mg2+ binding to specific sites, coupled with counter-ion release of monovalent cations, stimulate the formation of secondary and tertiary structures, leading to diverse pathways that include direct rapid folding and trapping in misfolded structures. In some molecules, key tertiary structural elements form when Mg2+ ions bind to specific RNA sites at the earliest stages of the folding, leading to specific collapse and rapid folding. In others, the formation of non-native base pairs, whose rearrangement is needed to reach the folded state, is the rate-limiting step. Escape from energetic traps, driven by thermal fluctuations, occurs readily. In contrast, the transition to the native state from long-lived topologically trapped native-like metastable states is extremely slow. Specific collapse and formation of energetically or topologically frustrated states occur early in the assembly process.

GHOST-NOT and GHOST-YES: Two programs for generating high-speed biosensors with randomized oligonucleotide binding sites with NOT or YES Boolean logic functions based on experimentally validated algorithms.

High-speed allosteric hammerhead ribozymes can be engineered to distinguish well between a perfectly matching effector and the nucleic acid sequences with a few mismatches under physiologically relevant conditions. Such ribozymes can be designed to control the expression of exogenous mRNAs and can be used to develop new gene therapies, including anticancer treatments. The in vivo selection of such ribozymes is a complicated and lengthy procedure with no guarantee of success. Thus, in silico selection of high-speed ribozymes can be employed using secondary RNA structure computation based on the partition function of the RNA folding in combination with random search algorithms. This approach has already been proven very accurate in designing allosteric hammerhead ribozymes. Herein, we present two programs for the computational design of allosteric ribozymes sensing randomized oligonucleotides based on the extended version of the hammerhead ribozyme. A Generator for High-speed Oligonucleotide Sensing allosteric ribozymes with NOT logic function (GHOST-NOT) and a Generator for High-speed Oligonucleotide Sensing allosteric ribozymes with YES logic function (GHOST-YES) for computational design of high-speed allosteric ribozymes are described. The allosteric hammerhead ribozymes had a high self-cleavage rate of about 1.8 per minute and were very selective in sensing an effector sequence.

Group I Intron as a Potential Target for Antifungal Compounds: Development of a Trans-Splicing High-Throughput Screening Strategy.

The search for safe and efficient new antifungal compounds for agriculture has led to more efforts in finding new modes of action. This involves the discovery of new molecular targets, including coding and non-coding RNA. Rarely found in plants and animals but present in fungi, group I introns are of interest as their complex tertiary structure may allow selective targeting using small molecules. In this work, we demonstrate that group I introns present in phytopathogenic fungi have a self-splicing activity in vitro that can be adapted in a high-throughput screening to find new antifungal compounds. Ten candidate introns from different filamentous fungi were tested and one group ID intron found in F. oxysporum showed high self-splicing efficiency in vitro. We designed the Fusarium intron to act as a trans-acting ribozyme and used a fluorescence-based reporter system to monitor its real time splicing activity. Together, these results are opening the way to study the druggability of such introns in crop pathogen and potentially discover small molecules selectively targeting group I introns in future high-throughput screenings.

Ribozyme-Catalyzed Late-Stage Functionalization and Fluorogenic Labeling of RNA.

Site-specific introduction of bioorthogonal handles into RNAs is in high demand for decorating RNAs with fluorophores, affinity labels or other modifications. Aldehydes represent attractive functional groups for post-synthetic bioconjugation reactions. Here, we report a ribozyme-based method for the synthesis of aldehyde-functionalized RNA by directly converting a purine nucleobase. Using the methyltransferase ribozyme MTR1 as an alkyltransferase, the reaction is initiated by site-specific N1 benzylation of purine, followed by nucleophilic ring opening and spontaneous hydrolysis under mild conditions to yield a 5-amino-4-formylimidazole residue in good yields. The modified nucleotide is accessible to aldehyde-reactive probes, as demonstrated by the conjugation of biotin or fluorescent dyes to short synthetic RNAs and tRNA transcripts. Upon fluorogenic condensation with a 2,3,3-trimethylindole, a novel hemicyanine chromophore was generated directly on the RNA. This work expands the MTR1 ribozyme's area of application from a methyltransferase to a tool for site-specific late-stage functionalization of RNA.

A ribozyme that uses lanthanides as cofactor.

To explore how an early, RNA-based life form could have functioned, in vitro selection experiments have been used to develop catalytic RNAs (ribozymes) with relevant functions. We previously identified ribozymes that use the prebiotically plausible energy source cyclic trimetaphosphate (cTmp) to convert their 5'-hydroxyl group to a 5'-triphosphate. While these ribozymes were developed in the presence of Mg2+, we tested here whether lanthanides could also serve as catalytic cofactors because lanthanides are ideal catalytic cations for this reaction. After an in vitro selection in the presence of Yb3+, several active sequences were isolated, and the most active RNA was analyzed in more detail. This ribozyme required lanthanides for activity, with highest activity at a 10:1 molar ratio of cTmp : Yb3+. Only the four heaviest lanthanides gave detectable signals, indicating a high sensitivity of ribozyme catalysis to the lanthanide ion radius. Potassium and Magnesium did not facilitate catalysis alone but they increased the lanthanide-mediated kOBS by at least 100-fold, with both K+ and Mg2+ modulating the ribozyme's secondary structure. Together, these findings show that RNA is able to use the unique properties of lanthanides as catalytic cofactor. The results are discussed in the context of early life forms.

RNA-Catalyzed RNA Ligation within Prebiotically Plausible Model Protocells.

Demonstrating RNA catalysis within prebiotically relevant models of primordial cells (protocells) remains a challenge in origins of life research. Fatty acid vesicles encapsulating genomic and catalytic RNAs (ribozymes) are attractive models for protocells; however, RNA catalysis has largely been incompatible with fatty acid vesicles due to their instability in the presence of Mg2+ at the concentrations required for ribozyme function. Here, we report a ribozyme that catalyzes template-directed RNA ligation at low Mg2+ concentrations and thus remains active within stable vesicles. Ribose and adenine, both prebiotically relevant molecules, were found to greatly reduce Mg2+ -induced RNA leakage from vesicles. When we co-encapsulated the ribozyme, substrate, and template within fatty acid vesicles, we observed efficient RNA-catalyzed RNA ligation upon subsequent addition of Mg2+ . Our work shows that RNA-catalyzed RNA assembly can occur efficiently within prebiotically plausible fatty acid vesicles and represents a step toward the replication of primordial genomes within self-replicating protocells.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology.

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4 kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.